Image-based multiplex immune profiling of cancer tissues: translational implications. A report of the International Immuno-oncology Biomarker Working Group on Breast Cancer.

Recent advances in the field of immuno-oncology have brought transformative changes in the management of cancer patients. The immune profile of tumours has been found to have key value in predicting disease prognosis and treatment response in various cancers. Multiplex immunohistochemistry and immunofluorescence have emerged as potent tools for the simultaneous detection of multiple protein biomarkers in a single tissue section, thereby expanding opportunities for molecular and immune profiling while preserving tissue samples. By establishing the phenotype of individual tumour cells when distributed within a mixed cell population, the identification of clinically relevant biomarkers with high-throughput multiplex immunophenotyping of tumour samples has great potential to guide appropriate treatment choices. Moreover, the emergence of novel multi-marker imaging approaches can now provide unprecedented insights into the tumour microenvironment, including the potential interplay between various cell types. However, there are significant challenges to widespread integration of these technologies in daily research and clinical practice. This review addresses the challenges and potential solutions within a structured framework of action from a regulatory and clinical trial perspective. New developments within the field of immunophenotyping using multiplexed tissue imaging platforms and associated digital pathology are also described, with a specific focus on translational implications across different subtypes of cancer. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Clinical Impact of New Treatment Strategies for HER2-Positive Metastatic Breast Cancer Patients with Resistance to Classical Anti-HER Therapies.

HER2-positive breast cancer accounts for 15-20% of all breast cancer cases. This subtype is characterized by an aggressive behavior and poor prognosis. Anti-HER2 therapies have considerably improved the natural course of the disease. Despite this, relapse still occurs in around 20% of patients due to primary or acquired treatment resistance, and metastasis remains an incurable disease. This article reviews the main mechanisms underlying resistance to anti-HER2 treatments, focusing on newer HER2-targeted therapies. The progress in anti-HER2 drugs includes the development of novel antibody-drug conjugates with improvements in the conjugation process and novel linkers and payloads. Moreover, trastuzumab deruxtecan has enhanced the efficacy of trastuzumab emtansine, and the new drug trastuzumab duocarmazine is currently undergoing clinical trials to assess its effect. The combination of anti-HER2 agents with other drugs is also being evaluated. The addition of immunotherapy checkpoint inhibitors shows some benefit in a subset of patients, indicating the need for useful biomarkers to properly stratify patients. Besides, CDK4/6 and tyrosine kinase inhibitors are also included in the design of new treatment strategies. Lapitinib, neratinib and tucatinib have been approved for HER2-positive metastasis patients, however clinical trials are currently ongoing to optimize combined strategies, to reduce toxicity, and to better define the useful setting. Clinical research should be strengthened along with the discovery and validation of new biomarkers, as well as a deeper understanding of drug resistance and action mechanisms.

Spatial analyses of immune cell infiltration in cancer: current methods and future directions: A report of the International Immuno-Oncology Biomarker Working Group on Breast Cancer.

Modern histologic imaging platforms coupled with machine learning methods have provided new opportunities to map the spatial distribution of immune cells in the tumor microenvironment. However, there exists no standardized method for describing or analyzing spatial immune cell data, and most reported spatial analyses are rudimentary. In this review, we provide an overview of two approaches for reporting and analyzing spatial data (raster versus vector-based). We then provide a compendium of spatial immune cell metrics that have been reported in the literature, summarizing prognostic associations in the context of a variety of cancers. We conclude by discussing two well-described clinical biomarkers, the breast cancer stromal tumor infiltrating lymphocytes score and the colon cancer Immunoscore, and describe investigative opportunities to improve clinical utility of these spatial biomarkers. © 2023 The Pathological Society of Great Britain and Ireland.

Pitfalls in machine learning-based assessment of tumor-infiltrating lymphocytes in breast cancer: A report of the International Immuno-Oncology Biomarker Working Group on Breast Cancer

The clinical significance of the tumor-immune interaction in breast cancer is now established, and tumor-infiltrating lymphocytes (TILs) have emerged as predictive and prognostic biomarkers for patients with triple-negative (estrogen receptor, progesterone receptor, and HER2-negative) breast cancer and HER2-positive breast cancer. How computational assessments of TILs might complement manual TIL assessment in trial and daily practices is currently debated. Recent efforts to use machine learning (ML) to automatically evaluate TILs have shown promising results. We review state-of-the-art approaches and identify pitfalls and challenges of automated TIL evaluation by studying the root cause of ML discordances in comparison to manual TIL quantification. We categorize our findings into four main topics: (1) technical slide issues, (2) ML and image analysis aspects, (3) data challenges, and (4) validation issues. The main reason for discordant assessments is the inclusion of false-positive areas or cells identified by performance on certain tissue patterns or design choices in the computational implementation. To aid the adoption of ML for TIL assessment, we provide an in-depth discussion of ML and image analysis, including validation issues that need to be considered before reliable computational reporting of TILs can be incorporated into the trial and routine clinical management of patients with triple-negative breast cancer. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

miR-146a-5p Promotes Angiogenesis and Confers Trastuzumab Resistance in HER2+ Breast Cancer.

Trastuzumab treatment has significantly improved the prognosis of HER2-positive breast cancer patients. Despite this, resistance to therapy still remains the main clinical challenge. In order to evaluate the implication of microRNAs in the trastuzumab response, we performed a microRNA array in parental and acquired trastuzumab-resistant HER2-positive breast cancer cell lines. Our results identified miR-146a-5p as the main dysregulated microRNA. Interestingly, high miR-146a-5p expression in primary tumor tissue significantly correlated with shorter disease-free survival in HER2-positive breast cancer patients. The gain- and loss-of-function of miR-146a-5p modulated the response to trastuzumab. Furthermore, the overexpression of miR-146a-5p increased migration and angiogenesis, and promoted cell cycle progression by reducing CDKN1A expression. Exosomes from trastuzumab-resistant cells showed a high level of miR-146a-5p expression compared with the parental cells. In addition, the co-culture with resistant cells' exosomes was able to decrease in sensitivity and increase the migration capacities in trastuzumab-sensitive cells, as well as angiogenesis in HUVEC-2 cells. Collectively, these data support the role of miR-146a-5p in resistance to trastuzumab, and demonstrate that it can be transferred by exosomes conferring resistance properties to other cells.

High VEGFR3 Expression Reduces Doxorubicin Efficacy in Triple-Negative Breast Cancer.

Due to the lack of specific targets, cytotoxic chemotherapy still represents the common standard treatment for triple-negative breast patients. Despite the harmful effect of chemotherapy on tumor cells, there is evidence that treatment could modulate the tumor microenvironment in a way favoring the propagation of the tumor. In addition, the lymphangiogenesis process and its factors could be involved in this counter-therapeutic event. In our study, we have evaluated the expression of the main lymphangiogenic receptor VEGFR3 in two triple-negative breast cancer in vitro models, resistant or not to doxorubicin treatment. The expression of the receptor, at mRNA and protein levels, was higher in doxorubicin-resistant cells than in parental cells. In addition, we confirmed the upregulation of VEGFR3 levels after a short treatment with doxorubicin. Furthermore, VEGFR3 silencing reduced cell proliferation and migration capacities in both cell lines. Interestingly, high VEGFR3 expression was significantly positively correlated with worse survival in patients treated with chemotherapy. Furthermore, we have found that patients with high expression of VEGFR3 present shorter relapse-free survival than patients with low levels of the receptor. In conclusion, elevated VEGFR3 levels correlate with poor survival in patients and with reduced doxorubicin treatment efficacy in vitro. Our results suggest that the levels of this receptor could be a potential marker of meager doxorubicin response. Consequently, our results suggest that the combination of chemotherapy and VEGFR3 blockage could be a potentially useful therapeutic strategy for the treatment of triple-negative breast cancer.

Role of SALL4 in HER2+ Breast Cancer Progression: Regulating PI3K/AKT Pathway.

Treatment for the HER2+ breast cancer subtype is still unsatisfactory, despite breakthroughs in research. The discovery of various new molecular mechanisms of transcription factors may help to make treatment regimens more effective. The transcription factor SALL4 has been related to aggressiveness and resistance therapy in cancer. Its molecular mechanisms and involvement in various signaling pathways are unknown in the HER2+ breast cancer subtype. In this study, we have evaluated the implication of SALL4 in the HER2+ subtype through its expression in patients' samples and gain and loss of function in HER2+ cell lines. We found higher SALL4 expression in breast cancer tissues compared to healthy tissue. Interestingly, high SALL4 expression was associated with disease relapse and poor patient survival. In HER2+ cell lines, transient overexpression of SALL4 modulates PI3K/AKT signaling through regulating PTEN expression and BCL2, which increases cell survival and proliferation while reducing the efficacy of trastuzumab. SALL4 has also been observed to regulate the epithelial-mesenchymal transition and stemness features. SALL4 overexpression significantly reduced the epithelial markers E-cadherin, while it increased the mesenchymal markers ß-catenin, vimentin and fibronectin. Furthermore, it has been also observed an increased expression of MYC, an essential transcription factor for regulating epithelial-mesenchymal transition and/or cancer stem cells. Our study demonstrates, for the first time, the importance of SALL4 in the HER2+ subtype and partial regulation of trastuzumab sensitivity. It provides a viable molecular mechanism-driven therapeutic strategy for an important subset of HER2-overexpressing patients whose malignancies are mediated by SALL4 expression.

Targeted immunotherapy against distinct cancer-associated fibroblasts overcomes treatment resistance in refractory HER2+ breast tumors.

About 50% of human epidermal growth factor receptor 2 (HER2)+ breast cancer patients do not benefit from HER2-targeted therapy and almost 20% of them relapse after treatment. Here, we conduct a detailed analysis of two independent cohorts of HER2+ breast cancer patients treated with trastuzumab to elucidate the mechanisms of resistance to anti-HER2 monoclonal antibodies. In addition, we develop a fully humanized immunocompetent model of HER2+ breast cancer recapitulating ex vivo the biological processes that associate with patients' response to treatment. Thanks to these two approaches, we uncover a population of TGF-beta-activated cancer-associated fibroblasts (CAF) specific from tumors resistant to therapy. The presence of this cellular subset related to previously described myofibroblastic (CAF-S1) and podoplanin+ CAF subtypes in breast cancer associates with low IL2 activity. Correspondingly, we find that stroma-targeted stimulation of IL2 pathway in unresponsive tumors restores trastuzumab anti-cancer efficiency. Overall, our study underscores the therapeutic potential of exploiting the tumor microenvironment to identify and overcome mechanisms of resistance to anti-cancer treatment.

Targeting HER2-AXL heterodimerization to overcome resistance to HER2 blockade in breast cancer.

Anti-HER2 therapies have markedly improved prognosis of HER2-positive breast cancer. However, different mechanisms play a role in treatment resistance. Here, we identified AXL overexpression as an essential mechanism of trastuzumab resistance. AXL orchestrates epithelial-to-mesenchymal transition and heterodimerizes with HER2, leading to activation of PI3K/AKT and MAPK pathways in a ligand-independent manner. Genetic depletion and pharmacological inhibition of AXL restored trastuzumab response in vitro and in vivo. AXL inhibitor plus trastuzumab achieved complete regression in trastuzumab-resistant patient-derived xenograft models. Moreover, AXL expression in HER2-positive primary tumors was able to predict prognosis. Data from the PAMELA trial showed a change in AXL expression during neoadjuvant dual HER2 blockade, supporting its role in resistance. Therefore, our study highlights the importance of targeting AXL in combination with anti-HER2 drugs across HER2-amplified breast cancer patients with high AXL expression. Furthermore, it unveils the potential value of AXL as a druggable prognostic biomarker in HER2-positive breast cancer.

Histological Grade and Tumor Stage Are Correlated with Expression of Receptor Activator of Nuclear Factor Kappa b (Rank) in Epithelial Ovarian Cancers.

The receptor activator of nuclear factor kappa B (RANK) is becoming recognized as a master regulator of tumorigenesis, yet its role in gynecological cancers remains mostly unexplored. We investigated whether there is a gradation of RANK protein and mRNA expression in epithelial ovarian cancer (EOC) according to malignancy and tumor staging. Immunohistochemical expression of RANK was examined in a cohort of 135 (benign n = 29, borderline n= 23 and malignant n = 83) EOCs. Wild type and truncated RANK mRNA isoform quantification was performed in a cohort of 168 (benign n = 26, borderline n = 13 and malignant n = 129) EOCs. RANK protein and mRNA values were increased in malignant vs. benign or borderline conditions across serous, mucinous and endometrioid cancer subtypes. Additionally, a trend of increased RANK values with staging was observed for the mucinous and serous histotype. Thus, increased expression of RANK appears associated with the evolution of disease to the onset of malignancy in EOC. Moreover, in some EOC histotypes, RANK expression is additionally associated with clinicopathological markers of tumor aggressiveness, suggesting a role in further progression of tumor activity.

Is Axillary Assessment of Ductal Carcinoma In Situ of the Breast Necessary in All Cases?

BACKGROUND: Staging of the axilla in women with ductal carcinoma in situ (DCIS) is a point of controversy. We aimed to assess whether there is a group of patients in whom axillary assessment can be avoided and whether the likelihood of underdiagnosis of infiltrating carcinoma is sufficient to justify this evaluation. MATERIALS AND METHODS: This was a multicenter, prospective, observational study of patients who were operated on between 2008 and 2018 in three Spanish hospitals, with a diagnosis by radiological or excisional biopsy of DCIS and clinically and radiologically negative axilla. RESULTS: A total of 530 patients with a preoperative diagnosis of DCIS were studied. An axillary assessment was performed in 77% of the patients. In 397 patients, selective sentinel lymph node biopsy was performed. Axillary involvement was found in 7.2% of all patients, which dropped to 2.15% if we only included DCIS diagnosed after a definitive anatomical pathology analysis. Underdiagnosis was correlated with the type of biopsy performed: the risk was 1.34 times as high if the biopsy was performed with a core needle. The risk of lymph node metastasis was higher when there was lymphovascular invasion and when mastectomy was performed. CONCLUSIONS: We propose an axilla management algorithm in patients with a preoperative diagnosis of DCIS. The patients who would benefit from sentinel lymph node biopsy would be those who are not candidates for breast-conserving surgery, those with a BIRADS 5 lesion biopsied by core-needle biopsy, and those whose definitive diagnosis is lymphovascular invasion.

First Nationwide Molecular Screening Program in Spain for Patients With Advanced Breast Cancer: Results From the AGATA SOLTI-1301 Study.

BACKGROUND: The SOLTI-1301 AGATA study aimed to assess the feasibility of a multi-institutional molecular screening program to better characterize the genomic landscape of advanced breast cancer (ABC) and to facilitate patient access to matched-targeted therapies in Spain. METHODS: DNA sequencing of 74 cancer-related genes was performed using FFPE tumor samples in three different laboratories with three different gene panels. A multidisciplinary advisory board prospectively recommended potential targeted treatments. The primary objective was to determine the success of matching somatic DNA alteration to an experimental drug/drug class. RESULTS: Between September 2014 and July 2017, 305 patients with ABC from 10 institutions were enrolled. Tumor sequencing was successful in 260 (85.3%) patients. Median age was 54 (29-80); most tumors were hormone receptor-positive/HER2-negative (74%), followed by triple-negative (14.5%) and HER2-positive (11.5%). Ninety-seven (37%) tumor samples analyzed proceeded from metastatic sites. Somatic mutations were identified in 163 (62.7%) patients, mostly in PIK3CA (34%), TP53 (22%), AKT1 (5%), ESR1 (3%), and ERBB2 (3%) genes. Significant enrichment of AKT1 mutation was observed in metastatic versus primary samples (9% vs. 2%; p=0.01). Genome-driven cancer therapy was recommended in 45% (n=116) of successfully screened patients, 11% (n=13) of whom finally received it. Among these patients, 46.2% had a PFS of ≥6 months on matched therapy. CONCLUSIONS: AGATA is the first nationwide molecular screening program carried out in Spain and we proved that implementing molecular data in the management of ABC is feasible. Although these results are promising, only 11% of the patients with genome-driven cancer therapy received it.

Interobserver variability in the assessment of stromal tumor-infiltrating lymphocytes (sTILs) in triple-negative invasive breast carcinoma influences the association with pathological complete response: the IVITA study.

High stromal tumor-infiltrating lymphocytes (sTILs) in triple-negative breast cancer (TNBC) are associated with pathological complete response (pCR) after neoadjuvant chemotherapy (NAC). Histopathological assessment of sTILs in TNBC biopsies is characterized by substantial interobserver variability, but it is unknown whether this affects its association with pCR. Here, we aimed to investigate the degree of interobserver variability in an international study, and its impact on the relationship between sTILs and pCR. Forty pathologists assessed sTILs as a percentage in digitalized biopsy slides, originating from 41 TNBC patients who were treated with NAC followed by surgery. Pathological response was quantified by the MD Anderson Residual Cancer Burden (RCB) score. Intraclass correlation coefficients (ICCs) were calculated per pathologist duo and Bland-Altman plots were constructed. The relation between sTILs and pCR or RCB class was investigated. The ICCs ranged from -0.376 to 0.947 (mean: 0.659), indicating substantial interobserver variability. Nevertheless, high sTILs scores were significantly associated with pCR for 36 participants (90%), and with RCB class for eight participants (20%). Post hoc sTILs cutoffs at 20% and 40% resulted in variable associations with pCR. The sTILs in TNBC with RCB-II and RCB-III were intermediate to those of RCB-0 and RCB-I, with lowest sTILs observed in RCB-I. However, the limited number of RCB-I cases precludes any definite conclusions due to lack of power, and this observation therefore requires further investigation. In conclusion, sTILs are a robust marker for pCR at the group level. However, if sTILs are to be used to guide the NAC scheme for individual patients, the observed interobserver variability might substantially affect the chance of obtaining a pCR. Future studies should determine the 'ideal' sTILs threshold, and attempt to fine-tune the patient selection for sTILs-based de-escalation of NAC regimens. At present, there is insufficient evidence for robust and reproducible sTILs-guided therapeutic decisions.

Identification of a Two-MicroRNA Signature in Plasma as a Novel Biomarker for Very Early Diagnosis of Breast Cancer.

The early diagnosis of breast cancer is essential to improve patients' survival rate. In this context, microRNAs have been described as potential diagnostic biomarkers for breast cancer. Particularly, circulating microRNAs have a strong value as non-invasive biomarkers. Herein, we assessed the potential of a microRNA signature based on miR-30b-5p and miR-99a-5p levels in plasma as a diagnostic biomarker for breast cancer. This two-microRNA signature was constructed by Principal Component Analysis and its prognostic value was assessed in a discovery cohort and blindly validated in a second cohort from an independent institution. ROC curve analysis and biomarker performance parameter evaluation demonstrated that our proposed signature presents a high value as a non-invasive biomarker for very early detection of breast cancer. In addition, pathway enrichment analysis identified three of the well-known pathways involved in cancer as targets of the two microRNAs.

Sentinel lymph node BIOPSY after neoadjuvant therapy in breast cancer patients with lymph node involvement at diagnosis. Could wire localization of clipped node improve our results?

INTRODUCTION: Sentinel lymph node biopsy (SLNB) after neoadjuvant therapy (NAT) in node-positive (N+) breast cancer patients at diagnosis remains a controversial issue, with no consensus on implementation or safety. OBJECTIVES: We sought to assess the accuracy of SLNB after NAT in biopsy-proven N+ cases at diagnosis and the efficacy and accuracy of wire localization of the clipped node to improve results. MATERIAL AND METHODS: A cross-sectional diagnostic technique validation study in N+ patients following NAT was performed. The biopsy-proven affected lymph node was clipped at diagnosis. SLNB and axillary lymph node dissection (ALND) were performed in cases of clinical-radiological lymph node response after NAT. For the purposes of our study we added wire localization of the clipped node. RESULTS: 103 patients were included (mean age, 54.4 years [± 12.7]). Wire marking was performed in 28 cases. The overall identification rate (IR) of SLN was 81.6%. The median number of nodes removed was 2 (range 2). The overall false negative rate (FNR) was 6.1%. Sensitivity and overall accuracy were 93.9% and 95.2%, respectively (area under curve 0.97). In the double-marked (clip and wire) group the FNR decreased to 0% and accuracy was 100%. Axillary pathologic complete response was observed in 24.3% of cases. CONCLUSIONS: SLNB is useful in node-positive patients at diagnosis who respond to NAT. Combining this with preoperative wire localization of the biopsied lymph node reduces the FNR without increasing the number of complications.

Aurora kinases in ovarian cancer.

Aurora kinases (AURK) are key regulators of the mitotic spindle formation. AURK is frequently overexpressed in ovarian cancer and this overexpression has been frequently associated with prognosis in these tumours. Interestingly, AURK have been shown to interact with DNA repair mechanisms and other cell cycle regulators. These functions have brought light to Aurora family as a potential target for anticancer therapy. In the last years, two clinical trials with different AURK inhibitors have shown activity in epithelial and clear-cell ovarian cancer. Although there is a lack of predictive factors of AURK inhibition activity, recent trials have identified some candidates. This review will focus in the functions of the AURK family, its role as prognostic factor in epithelial ovarian cancer and potential clinical implications.

MicroRNA-33b Suppresses Epithelial-Mesenchymal Transition Repressing the MYC-EZH2 Pathway in HER2+ Breast Carcinoma.

Downregulation of miR-33b has been documented in many types of cancers and is being involved in proliferation, migration, and epithelial-mesenchymal transition (EMT). Furthermore, the enhancer of zeste homolog 2-gene (EZH2) is a master regulator of controlling the stem cell differentiation and the cell proliferation processes. We aim to evaluate the implication of miR-33b in the EMT pathway in HER2+ breast cancer (BC) and to analyze the role of EZH2 in this process as well as the interaction between them. miR-33b is downregulated in HER2+ BC cells vs healthy controls, where EZH2 has an opposite expression in vitro and in patients' samples. The upregulation of miR-33b suppressed proliferation, induced apoptosis, reduced invasion, migration and regulated EMT by an increase of E-cadherin and a decrease of ß-catenin and vimentin. The silencing of EZH2 mimicked the impact of miR-33b overexpression. Furthermore, the inhibition of miR-33b induces cell proliferation, invasion, migration, EMT, and EZH2 expression in non-tumorigenic cells. Importantly, the Kaplan-Meier analysis showed a significant association between high miR-33b expression and better overall survival. These results suggest miR-33b as a suppressive miRNA that could inhibit tumor metastasis and invasion in HER2+ BC partly by impeding EMT through the repression of the MYC-EZH2 loop.

The GATA3 X308_Splice breast cancer mutation is a hormone context-dependent oncogenic driver.

As the catalog of oncogenic driver mutations is expanding, it becomes clear that alterations in a given gene might have different functions and should not be lumped into one class. The transcription factor GATA3 is a paradigm of this. We investigated the functions of the most common GATA3 mutation (X308_Splice) and five additional mutations, which converge into a neoprotein that we called "neoGATA3," associated with excellent prognosis in patients. Analysis of available molecular data from >3000 breast cancer patients revealed a dysregulation of the ER-dependent transcriptional response in tumors carrying neoGATA3-generating mutations. Mechanistic studies in vitro showed that neoGATA3 interferes with the transcriptional programs controlled by estrogen and progesterone receptors, without fully abrogating them. ChIP-Seq analysis indicated that ER binding is reduced in neoGATA3-expressing cells, especially at distal regions, suggesting that neoGATA3 interferes with the fine tuning of ER-dependent gene expression. This has opposite outputs in distinct hormonal context, having pro- or anti-proliferative effects, depending on the estrogen/progesterone ratio. Our data call for functional analyses of putative cancer drivers to guide clinical application.

FGFR4 regulates tumor subtype differentiation in luminal breast cancer and metastatic disease.

Mechanisms driving tumor progression from less aggressive subtypes to more aggressive states represent key targets for therapy. We identified a subset of luminal A primary breast tumors that give rise to HER2-enriched (HER2E) subtype metastases, but remain clinically HER2 negative (cHER2-). By testing the unique genetic and transcriptomic features of these cases, we developed the hypothesis that FGFR4 likely participates in this subtype switching. To evaluate this, we developed 2 FGFR4 genomic signatures using a patient-derived xenograft (PDX) model treated with an FGFR4 inhibitor, which inhibited PDX growth in vivo. Bulk tumor gene expression analysis and single-cell RNA sequencing demonstrated that the inhibition of FGFR4 signaling caused molecular switching. In the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) breast cancer cohort, FGFR4-induced and FGFR4-repressed signatures each predicted overall survival. Additionally, the FGFR4-induced signature was an independent prognostic factor beyond subtype and stage. Supervised analysis of 77 primary tumors with paired metastases revealed that the FGFR4-induced signature was significantly higher in luminal/ER+ tumor metastases compared with their primaries. Finally, multivariate analysis demonstrated that the FGFR4-induced signature also predicted site-specific metastasis for lung, liver, and brain, but not for bone or lymph nodes. These data identify a link between FGFR4-regulated genes and metastasis, suggesting treatment options for FGFR4-positive patients, whose high expression is not caused by mutation or amplification.